Velocity anisotropy and the antiferromagnetic instability in the QED% theory of 

underdoped cuprates 
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We consider the effective 2+1 dimensional electrodynamics (QED3) of low-energy quasiparti- 
cles coupled to fluctuating vortex loops in the d-wave superconductor, with the velocity (Dirac) 
anisotropy vf/va 7^ 1. Fhis theory should be relevant to the quantum superconductor- insulator 
transition in underdoped cuprates. Working in the customary large-iV approximation, we find that 
weak anisotropy is a marginally irrelevant perturbation to the Lorentz invariant QED3, and that 
the critical number N c of Dirac fields below which the theory suffers the antiferromagnetic (chiral) 
instability stays the same. 



Two states of matter ubiquitous in all high- 
temperature superconductors are the insulating antifer- 
romagnet and the d-wave superconductor. Recently one 
of us has shown Q that, while experience has taught 
us that to derive superconductivity from antiferromag- 
netism is rather difficult, the opposite appears possible: 
d-wave superconductors (dSC) at T = have an inherent 
instability towards antifcrromagnetism once they loose 
their phase coherence. The low-energy theory that de- 
scribes this quantum phase transition is the anisotropic 
2+1 dimensional quantum electrodynamics (QED3), in 
which the two (N = 2) Dirac fermions correspond to 
neutral spinons near the two pairs of diagonally opposed 
nodes in the d-wave state, and the massless gauge field Q 
describes the interaction with the unbound vortex loops. 
Such a theory in the isotropic limit is well known to 
have a " chiral" instability for the number of Dirac com- 
ponents N < N c , with N c w 3 ||. In the context of 
d-wave superconductivity, this instability corresponds to 
the dynamically generated staggered potential felt by the 
original electrons, i.e. a weak spin density-wave (SDW) 
order @. The SDW obtained this way is then enhanced 
by Coulomb repulsion, and is assumed to be smoothly 
connected to the commensurate antiferromagnet right at 
half filling. Another crucial result of this approach is 
that such spontaneously generated SDW automatically 
confines neutral spin-1/2 excitations (spinons) of the su- 
perconductor. If, on the other hand, the critical number 
of components were N c < 2, the ground state would be 
qualitatively different: the phase incoherent dSC would 
not show any kind of order at T — 0, and the ground 
state would be analogous to the "algebraic Fermi (spin) 
liquids" H] , {§ , invoked in literature to describe the pseu- 
dogap regime in cuprates. Evidently, it is of some impor- 
tance to determine in which way the perturbations to the 
Lorentz invariant QED3, relevant to cuprates, affect the 
critical number of components for the SDW instability. 

One such perturbation that is quite significant in high- 
temperature superconductors is the large observed differ- 
ence between the two characteristic velocities that appear 
in the linearized Hamiltonian for the d-wave quasiparti- 
cles, vf/va ~ 10 Q. Here vp is the Fermi velocity, and 
va is related to the amplitude of the superconducting or- 



der parameter. Similar velocity anisotropy also appears 
in the other incarnations of the QED3 in the theories of 
high-T c superconductivity Q, and is in fact generic to 
the effective Dirac theories describing condensed matter 
systems 0, which have no reason to be Lorentz invariant. 
Because of the coupling to the gauge-field such anisotropy 
can not be simply rescaled out, and is in fact precisely a 
marginal perturbation to the Lorentz invariant QED3 by 
simple power counting. In this paper we study the effect 
of weak anisotropy of this type on the critical number of 
components N c in the effective QED3. 

Working in a large-N approximation, we find two re- 
sults. Firstly, that weak anisotropy becomes irrelevant 
to the leading order in 1/N. Secondly, we argue that N c 
is the same in the weakly anisotropic case as it is in the 
isotropic case. 

We calculate both the renormalized anisotropy and 
A N c by performing expansions in A' = A — 1 , 5' = S — 1 
and 1/N, where A = Vf/va and S = ^/vfva/c (with c 
being the propagation velocity of the gauge field), around 
the self-consistent approximation scheme of ||. Essen- 
tially, when we come to calculate N C1 this becomes an 
expansion m powers of 1/N° around N°, the value of N c 
obtained by || in the isotropic limit, as well as in powers 
of A' and 5' . In such an expansion we find that to lead- 
ing order in 1 /N® N c is an increasing function of weak 
anisotropy. Since we find that weak anisotropy is irrel- 
evant to 0(1 /N), when we rewrite our result in terms 
of the renormalized anisotropics Xr and Sr (to leading 
order in l/N®, this corresponds to setting A = Xr and 
5 = S R ) we find that N c = iV c °. 

The effective low-energy theory for the d-wave quasi- 
particles (i. e. spinons) at T = coupled to the topo- 
logical defects (vortex loops J8| ) in the superconducting 
phase can be written as £ = L\ + £2 + £a with Q: 



N/2 



7o(<9 - ieA ) + — 7i(<3i - ieA x ) 



+SV\^ 2 (d2 - ieA 2 ) 
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and C 2 is the same as C\, except that 1/vA <-> VA. 
We rescaled the coordinates so that 8 = ^vfva/c and 
A = vp/va- The 7-matrices are 70 = (i/v / 2)o'i <8> 
71 = (i/^/2)a 2 ®03, and 72 = — (i/V2)cr 2 ®ai, and satisfy 
the Clifford algebra {7^, 7„} = — <5 M „. <? are the standard 
Pauli matrices. The gauge field A represents the singular 
superconducting phase fluctations induced by the vortex 
loops, and VP fields are the neutral spin-1/2 excitations 
one can define near the nodes in the dSC. The coupling 
constant e ("charge") is proportional to the dual order 
parameter that signals the appearance of infinitely large 
vortex loops, i.e. the destruction of the phase coherence 
in the dSC jy, 0. We have assumed N/2 identical copies 
of each type of spinon field; the value of physical interest 
being N=2. Weak short-ranged interaction terms quartic 
in 4" are irrelevant by power counting and omitted in Eq. 
1. The last term in Eq. 1 is the gauge-fixing term, and as 
customary we will work in the Landau gauge e = 0. Due 
to the 1 / VA — ► \f\ exchange in L 2 it will prove necessary 
to keep both 8 and A anisotropics in the action, once 
renormalization effects are included. In the limit 8 1, 
A — > 1 the model simply reduces to the much studied 
TV-flavour QED 3 . 

In the calculation of the polarization tensor 11^ ^, for 
5 7^ 1 and A ^ 1, we use the bare spinon propagators and 
work only to leading order in 1/iV, as shown in Fig. I. In 
Appendix A we give explicit forms of and 11^. The 
full gauge field propagator obeys the following symmetry 
properties: 

A),o(Po,Pi,P2) = A),o(Po,P2,Pi) 
Di,i(jpo,Pi,P2) = D 2>2 (po,p 2 ,Pi) 
Do,i(po,Pi,P2) = D 0i2 (p ,p 2 ,px) 
Di, 2 (p ,Pi,P2) = D h2 (p ,p 2 ,pi). (2) 

In what follows, for the sake of clarity, we shall omit the 
calculational details and just outline our method. The in- 
terested reader is referred Appendix B where some of the 
more important details are given. We start this discus- 
sion by considering the form the full spinon propagator 
at the ith node should take 

S hR (k) = (Ei(fc) + A^fch^))' 1 . (3) 

Here, the repeated index, as usual, implies a sum over 
indices. Ej are the dynamically generated mass-gaps. We 
refer to A^ as the wave-function renormalization in k^ 
direction, for the spinon at the ith node. Using our form 
for the full gauge-field propagator and the above form 
for the spinon propagator we then proceed to construct 
the approximate Schwinger-Dyson equations for the full 
spinon propagators as shown diagrammatically in Fig. 2. 

From Eq. 2 and the form of the Schwinger-Dyson equa- 
tions, it is easy to show that the following symmetry 
properties must also apply for the quantities we are in- 
terested in: 

Ei(feo) h,k 2 ) = E 2 (fc , k 2 ,ki) 
A lii (k ,ki,k 2 ) = A 2 j(k ,k 2 ,k 1 ). (4) 



We therefore need consider only type 1 spinons and here- 
after drop the spinon index. At this point it is also con- 
venient for us to introduce the rcnormalizcd anisotropics 
A;? = A 2 (0)/Ai(0) and Sr = y /A 1 (0)A 2 (0)/A (0) . 

From the Schwinger-Dyson equations we are able to 
generate four independent equations which we shall anal- 
yse. For the purposes of illustration, we will give here, 
only, the mass-gap equation. The other equations corre- 
spond to Aq, Ai and A2 from which we derive expressions 
for Aij and 8r; these are given in Appendix B. The mass 
equation takes the form 

where 

8 2 

U(k) = k 2 a + \8 2 k\ + —kf + EiOfc) 2 
A 

fC(k) = D 0i0 (k) + yDi,i(£) + 8 2 XD 2 , 2 (k) (6) 

with q = k — p, k = (ko,ki,k 2 ), and a — Ne 2 . These 
equations are difficult to use in their present form, due to 
the cumbersome form of K,. Therefore, we limit ourselves 
to the case of weak anisotropy, and treat 8' = 6 — 1 and 
small A' = A — I as small parameters. 

By expanding the mass-gap equation up to second 
order in both 8' and A' we may write N c in following 
form that depends on a rescaled self-energy E r (p/E) = 
E(p)/E, where E = E(0), 

AT C = / 7 0E r (fc)[F^(fc,d) + A'F 1 -°(fc,O) 

+8'F - 1 (k,6) + S'\'F 1 > 1 {k 1 d) 
+8' 2 F°- 2 {k, 0) + X /2 F°- 2 (k, 0)] (7) 

where the forms for F°- Q -F ' 2 are given in Appendix B, 
but with E(fc) replaced with E r (fc). 

Next we further approximate E r (fc) with its form in 
the isotropic limit, E°(fc), consistent to leading order in 
the l/iVj expansion which we employ. The 0(1/AT°) cor- 
rection to E r is then given by 

+S'F°' 1 (k, p; E°) + S'X'F^ik, p; E°) 
+S' 2 F°' 2 (k,p; E°) + 8' 2 F 2 >°(k,p; E°)]. (8) 

By replacing E°(fc) with E°(fc) + AE r (fc) in the equa- 
tion for N c and expanding to first order in AE r one is 
able to generate an 0(1/A°) term in the expansion of 
A7V C ; which is finite. However, to be fully consistent in 
calculating the 0(1/ N^) correction in such a calculation, 
together with the correction AE, corrections arising from 
the vertex and wave- function renormalization must be in- 
cluded. Indeed, it is well known that these change the 
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value of JV° even in the fully isotropic case ||. In our 
calculation, all the terms that contribute to 0(1/ N®) are 
shown in Fig 3. Terms that contribute to AN C , gener- 
ated from the vertex correction and wave-function renor- 
malizations will contain logarithmic singularities. As we 
shall see, S' R and A^ will also contain logarithmic singu- 
larities. So we assume we may absorb the singularities in 
AA C into the renormalized quantities S' R and X' R , which 
would provide us with a well defined expansion. 

We now calculate the leading order term in this expan- 
sion. It is easy to show that both the coefficients of A' and 
A'<5' vanish under angular integration. We are then left 
with the problem of determining the Ej?. One cannot, 
simply, set E°(fc) = 1, for one would obtain N® = oo, 
contrary to the values obtained by direct solution of the 
mass-gap equation in the isotropic limit. in fact 

must be a monotonicaly decreasing function of rescaled 
momentum. Based on the solution in jj| and the fact that 
that E°(0) = 1, we propose a simple ansatz for S°(fc) at 



Vk ' 
A 



(9) 



By inserting this form into the mass-gap equation Eq. 5 
in the isotropic limit one is able to find equations that 
fix A and determine N®. The points p = and p — A are 
only considered; at which we may set E°, to 1 and 1/2, 
respectively. The solution to these equations is then 

A = 5.60 and N° = 3.19. This is very close to 
N° = 3.24 f|, so the ansatz Eq. 9 will be sufficient 
for our purposes. We then use this to determine AA C as 
a function of bare anisotropics: 



AA C = 0.4295' + 0.5665' 2 + 0.406A' : 



(10) 



Furthermore, we are able to show that both the 8' and 
the A' 2 coefficients are positive definite, for all possible 
forms of Xy(fc). Although we are unable to prove this 
generally for the S' 2 term, we have been able to show that 
it remains positive over all chosen values of the parameter 
A; not just for A = 5.60. 

As argued above, we then take AA^ C to be given by the 
above expression, but with the bare anisotropies replaced 
with the renormalized values. Since the latter differ from 
the bare ones by terms of 0(1 /N), this is certainly cor- 
rect to the zeroth order in the expansion in 1/N. We 
are assuming, however, that by doing this we are effec- 
tively including the leading logarithmic divergences from 
the vertex corrections in Fig. 3., which are formally of 
the next order in 1/N. One justification for this assump- 
tion will be that by doing this we will find that N c be- 
comes independent of the parameter <5; as it should be, 
since we could have rescaled the theory so to place 6 
entirely into the gauge-field propagator and the gauge 
fixing term. Since the former is irrelevant compared to 
the fermion polarization, gauge invariance ensures any 
dependence on 5 in Eq. 1 does not effect physical quan- 
tities. The other justification is that one could expect 



a quantity like N c to be universal; i.e. independent of 
all weak irrelevant perturbations. In the example of the 
chiral instability in presence of weak quartic interactions, 
which are irrelevant N c , indeed, does not change [jlo| . It 
is also easy to check that in the isotropic QED3 including 
higher derivatives in the fermionic or gauge fields can not 
affect N°, to the leading order in 1/N. If, indeed, N c is 
universal we can fully replace the bare anisotropies with 
renormalized ones and show that N c will be independent 
of anisotropy, since we will ultimately find anisotropy to 
be an irrelevant perturbation. 

We therefore next determine X R = 1 + X' R and 5 R = 
1 + S' R from our approximate Schwinger-Dyson equations. 
For the equations for A, S' R and A^ we find forms similar 
to that of Eq. 8, when we expand in powers of 5' and 
A'. In the limit E/a — > 0, certain terms develop singu- 
larities, which dominate the behaviour of Ag, 6r and Xr. 
Expressed in terms of these singularities we find 



An = 1 



M 

N 



3vr 2 A 



1 




8' 



A'. (11) 



A, B, and C, are finite numerical constants. From these 
we may read off the flow of anisotropies to order 0(1/ N): 



ds 
ds 



64 
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5', 



5tt 2 A A ^' 



(12) 



Weak S' R and X' R are thus both 



where s = ln(a/E 
irrelevant as E — > and flow to zero. This agrees with 
the conclusion of [|1J . From this and the Eq. 9 we deduce 
that to leading order in 1/N°, AN C = 0. 

One should emphasise, for the above conclusion to hold 
to all orders in our expansion in N c it appears essential 
that N c is written in terms of renormalized anisotropies 
only, which would render it a universal number, depen- 
dent only on the infrared behaviour in the theory. Al- 
though one would expect this to be true in all cases, based 
on the examples given previously; it is by no means obvi- 
ous that our expectation is true in the case of anisotropy, 
or for that matter, other marginal perturbations. 

To summarize, we find that AN C = to leading order 
in 1/N®. This implies that weak velocity anisotropy does 
not effect N c and that one is free to set vf = «a = c = 1 
in the determination of all critical quantities. However, 
issues such as the universality of N c and the regime of 
strong anisotropy, need to be considered before one is 
certain that this is, indeed, the case beyond the limits 
of our calculation. We hope to address such problems in 
future work. 

This work is supported by NSERC of Canada and the 
Research Corporation. 
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APPENDIX A: FULL GAUGE FIELD 
PROPAGATOR 

To calculate the full gauge field propagator, we first 
need to calculate the polaziation function given in Fig. 
1. We first note that from the Feynman rules 

r d 3 k 

= £ J j^%,iSi,o(k)%,iSo(k + p) (Al) 

where 70,1 = 70, 71,1 = 6/VX71 , and 72,1 = <5VA~7 2 ; 
and also 7^,2 (A) = 7^,2(1/^). S^o is the bare spinon 
propagator for the ith node, and takes the form 



Si,o 



(A2) 



Here, as usual, a repeated index implies a sum over that 
index. It is easy to deduce from the calculation of II^„ 
in the isotropic limit what explicit form the elements of 



^ should take 



1 



n . = ni(p ,pi,P2) 

Il 1A (p) = Il 2i2 (p E ) = U 2 {p) 



'■ n 2 ,o(p*E) = popin 3 (p) = 

popi r 1 _a_ 

8 [APi P 2 

ni, 2 (p) = -pip 2 ^4(p) = 



P1P2S 



2 r 



1 



1 



(A3) 



Here, for convenience, we employ the notation: p = 

(P0,Pl,P2),PE = (j>0,P2,Pl), Pi = {pl + ^Pl + ^Xpl) 1 / 2 

and P 2 (A) = Pi(l/A). One may immeadiatly write down 
the expression for the full inverse gauge field propagator 

D£ = k 2 (V - (1 - e)^) + |H„ (A4) 

Straightforwardly although awkwardly, from Eq. A. 4 one 
is able to invert the inverse gauge field propagator, and 
so obtain expressions for the elements of the gauge-field 
propagator for arbitrarily covariant gauge. 



D ,o(P) 
Do,i(P) 



Di{p) D hl {p) = D 2 . 2 (pe) = D 2 (p) 
D Ape) = D 3 (p) D X AP) = Mp) 
Di{p) = V e (p)(D l AP) + zDiAp)) 



where 



Di.AP) 



2(2 

Pi [p 



Mpe)) + 



p\ (p 2 + |n 2 (p)) - aplp 2 2 Il 4 (p) 

d 1 Ap) = (pI+p\ + ^Mpe)) 



> \PI+Pl + ^(p) 



-Philip) -I 

p 2 (p 2 + ^U 2 (p E )) 



+ 



D2AP) 

P2 (p 2 + I n i(p)) ~ aplpl^iPE 

d 2 Ap) = (po + pi + fn 2 (p B )) 
x (pI+pI + f Si(#) -p 2 p 2 (|n 3 (p B ) - 1) 2 
d 3 ,Ap) = [plf (Ss(fe) + n 4 (p) - Ss(fe)) 
- (V + f n 2 (p £ ) 

r - / rv _ 

DzAp) =PoPi 



P0P1 



p- 2 [^ti 3 (p E )-i) (|n 4 (p)-i 



DaAp) = 



pI^{Mpe)+Mp)-Mp)) 



P1P2 



(p 2 + ^n l{ p)) 

^%{p E ) - l) 

(|n 4 (p) - 1) (pl+pl + 1 flxG?))] ,(A6) 



DaAp) = P0P1 Pa ( ^ n 3(p) - 1 



and 



where 



pl + e{p\ +pl+fli(p)) 
VAp) + eV e AP) + ^eAP) 



(A7) 



©e,l(^) = ^(p)A e (p E ) 

+A L (p £ )A e (p) - 2B L (p)B t (p) 
V e AP) = Mp)BM - BM 2 



(A8) 



and 



+P\{, 



ri[p 2 + ^2(p) 



p 2 + -U 1 (p)J -ap 2 p 2 U 3 (p) 
Bl{p) = P1P2 (p 2 + ^i(P)) 



+p%2 (Mp)-Mp)-Mpe)) 



A e (k) 



Pl+Pl+ 2 fl 2(fc) 



(pI+pI + 1 Si (£)) -pgp 2 (|ns(po,Pi,P2) - 1)' 
(A5) Be (^ = P1P2 [(|n 4 (p) - l) (p 2 + p 2 + |nr(p) 



-Po (fn 3 (p) - 1) (|n 3 (p £ ) 



(A9) 
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APPENDIX B: DETAILS IN DERIVING THE 
FORMULA FOR N c , 8 R AND \ R FROM 
SCHWINGER-DYSON EQUATIONS 

In what follows wc shall limit ourselves to the Lan- 
dau gauge. From Fig. 2 we may easily write down the 
expressions for the S-D equations 

d 3 k _ 



Again, it is easy to show from the symmetry properties 
of D^ v that ^1,1(^0,^1,^2) = A 2}i {k ,k 2l ki). For the 
moment we will concentrate on the mass-gap equation 
and calculate N c , before looking at wave- function renor- 
malization. 

Now, as discussed in the text, we perform an expansion 
in both 5' = S — 1 and A' = A — 1 which gives us up to 
second order 



Where Si,R and 5j,D arc 

s hR (k) = (Si(^) + AA^k)' 1 

S i>D {k) = (£*(£)+ TMJV)" 1 . (B2) 

By performing the trace of Eq. B.l we arrive at the 
mass-gap equations 



where 



+S'F°' 1 (k,p) + 8'X'F 1 - 1 (k,p) 
+S' 2 F°' 2 (k,p) + \' 2 F°> 2 (k,p)} (B8) 

F^(k,P) = -^ 

q 

16K_ 8Q_ 



^1,0 = 



F o,i = 



Pq q 3 I ' 
I6Q0 32if+ 



v " f d 3 k Si(fc) 



(B3) 



F i,i = 



Iq 3 Pq ' 

2K_ AK + K_\ 16 



J 2 



P 



1 



where 



8K-Q - 16g+Q_ 8q 2 Q_ - 8g 2 Q_ 

T" r 1 ) Q I 



LTi(fc) = fcg + A<5 2 fc 2 + fc 2 + Sl (£) 2 
A 

t/ 2 (fc) = k 2 + ^-k 2 + S 2 Xk 2 + E 2 (fc) 2 
A 

/Ci(fc) = A), (fc) + + * 2 ^2, 2 (fc) 



^2,0 = 

F 0,2 = 



J2 g 3 q sj 
,.2\ 1C f /-) /Q„2„2 o^)2 

7T J2 

AK\ K+\ 16 32^+Qo , 8g 2 Q 
3 



fej^ 16 if_Q_ (8g 2 g 2 - 2 Q 2 
g 7 2 <7 3 q b I 



P 



P 



Iq 3 



+ 



q b I 



^B9) 



where K+ = k\ + fcf, K_ = k 2 - Q+ = <? 2 + gf, 
MA) = °o,oW + <5 2 AL>i 4 (fc) + °—D 2 a{k). (B4) Q- = q\ - Qo = q 2 + q 2 , I(k) = k 2 + E(fc) 2 and 

<f = fc — p. By rescaling the self-energy E r (p') = E(p)/E, 
By using the symmetry properties given by Eq. 2, in the loop momenta k' = fc/E; and setting external momenta 
text, it is very easy to show that, indeed, Ei(fc , h,k 2 ) = P = we arrive at 7 g iven in the text • B Y performing 
E(fco,£>2,fci)- By multiplying Eq. B.l from the left with the expansion in 5' and A' to first order we find at an 
70 and then taking the trace we arrive at the equation expression for A 
for A lfi 



aMp) = 1 



a f d 3 k 2%(k,q) - fcoAC(g) 



A (p) = 1 



d 3 k 



(F°< (k,p)-k F '°(k,p)) 



where 



N Po J (2tt)3 



<5 2 fci 



, (B5) 



N Po J (2tt) 3 ' 

+X'(F^ (k,p)-k F 1 '°(k,p)) 

+S'(F^ 1 (k,p)-k F°' 1 (k,P)) 

(BIO) 



%{k,q) = k Q D t Aq) + ^ L D hl (q) + 5 2 Xk 2 D ia (q). (B6) where 

Similarly, by multiplying Eq. B.l by either 71 or 72 one 
obtains equations for Ai t i and Ai y2 , respectively 

All ffl = J-( 1 a f d 3 k 2T 1 (k,q)-k 1 IC(q)^ 



^'°(k,p) = 



(k Q+ - qo(qih + q 2 k 2 )) 16 



VA 



N Pl J (2tt) 3 



U(k) 



q 3 I{k) 

pi,0 _ ]QK- (k Q+ - qojqih + q 2 k 2 )) 
~ Pq 3 

1 16q (k 2 q 2 - fcigi) 
/ 3 



^0 



^oQ+ + {qih - q2k 2 )Q-q 



6 



2(fcigi + fc 2g2 ) go A 16 
q 3 J I 

32K + (k Q Q + - gp(gifci + g2 fc 2 )) 
7 2 9 3 



(Bll) 



For A\ we find 

+A'(F 1 1 ^°(fc,p)-fc lJ F 1 ^(fc,p)) 
+6'(F^ 1 (k,p)-k 1 F^(k,p))) A 



where 



(B12) 

pi,o 16JT- (fciQf - gi(gofco + g 2 fc 2 )) 
i ~ /2(? 3 

16 (fciQf + fc 2glg2 ) 
Iq 3 

k 2 qiq 2 Q+ _ fc Q_ go gi _ fci (Q+Q+ - 2 g g g 2 ) 



h2 (fc^o 2 - /, 



: 2 gig 2 



where 



(B14) 



f° °r£ -*> - ( fc 2Q+ - g 2 (gofco + gifci)) 16 
16if_ (fc 2 Q° 1 - g2 ( g0 fc + g 1 fc 1 )) 



f?' 1 = 



16 (fc 2 <ffi + fc lg2gl ) 

7 g 3 

fcigig2<3+ fcoQ-gog 2 fc 2 - 2 g p g |) 



2 (fc 2 g° 1 - fci g i g2 ) 



16 

T 



16 

T 

32K+ (hQf - gl ( go fc + g2 fc 2 )) 
p q 3 

(B13) 

where Q^ 2 = g Q + g 2 ■ For A 2 we find 

+A'(F 2 1 ' (fc,p)~fc 2 F 1 '°(fc,p)) 
+S'(F°'\k,p) - k 2 F°' 1 (k,p)j) Sy/X 



32K+ (k 2 Qf - q 2 (q k a + gl fcx)) 

/2g 3 

(B15) 

where Q+ 1 = g „ + g 2 . 

Now, we are interested in calculating the singular con- 
tributions to Xr and 8r in the limit £ — > 0. To do this 
we require the singular contributions to Aq(0), ^i(O) and 
yl 2 (0). In finding the logarithmically singular contribu- 
tions one may neglect the momentum dependence of £(p) 
and simply replace it with S(0). Then one is able to per- 
form the angular integrations. We find the following 



An M 1 - 



A 2 * 



3tt 2 N 
1 - 



37T 2 iV 



1 - 



1 



3vr 2 iV 



4 4 
7 ~ 5 
6Y 
5 
6V 
5 



5' 



In 



(I) 



ln(-))*VX. (B16) 



Where, we have neglected non-singular 0(1 /N) terms. 
Using the definitions for 5r and Xr given in the text, we 
thereby obtain Eq. 11 . 
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FIG. 3: Graphs that contribute to the 0(1/N°) correction to 
AN C . Closed spinon loops imply a sum over both nodes. 



FIG. 1: The polarisation tensor n Mi „ to the leading order in 
1/N. The thin straight lines with the numbers one and two 
represent bare spinon 1 and 2 propagators. 



FIG. 2: The Schwinger-Dyson equations. The full spinon 
propagators are represented by thick lines. Spinon propaga- 
tors that contain unrenormalized anisotropy (<5,A), but the full 
mass gap, are represented by double lines. The full gauge-field 
propagator is represented by a wavy line. 



